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A THEORY OF UNSTAGGERSD AIRFOIL CASCADES 
IN COMPRESSIBLE FLOW. 

By Robert A. Spurr a.n& H. Julian Allen 


SUMMARY 

By use of the methods of thin airfoil theory, which include 
effects of compressibility, relations are developed which permit 
the rapid determination of the pressure distribution over an 
unstaggered cascade of airfoils of a. given profile, and the deter- 
mination of the profile shape necessary to yield a given pressure 
distribution for small chord/gap rati os . For incompressible flow 
the results of the theory are compared with available examples 
obtained by the more exact method of conformal transformation. 
Although the theory is developed for small chord/gap ratios, these 
comparisons show that it may be extended to chord/gap ratios of 
order unity, at least for low— speed flows. Choking of cascades, 
a phenomenon of particular importance in compressor design, is 
considered. 


INTRODUCTION 

The wider use of gas turbines and other devices employing 
axial— flow compressors has increased the need for compressors with 
a high pressure rise per stage. In order to achieve this purpose, 

•it is necessary to use high velocity flows, thus increasing the 
possibility of losses through compression shock. A. method is there- 
fore desirable which will permit the design of compressor blades 
which have high critical compressibility speeds. This result can 
be accomplished if a cascade of airfoils representing the flow can 
be designed to give a desirable airfoil— section, pressure distribution. 

This report attacks a portion of the problem by finding the 
relation between the profile shape and pressure distribution over 
airfoils arranged in an unstaggered cascade through the use of the 
approximate methods of thin airfoil theory originally presented by 
Glauert in reference 1 and further developed by the NACA in 
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reference 2. The flow over an airfoil in cascade is related to that 
over a single airfoil in a free stream . The problem of finding the 
pressure distribution over an airfoil in cascade or the shape of 
an airfoil in cascade to give a required pressure distribution then 
reduces to the analogous problem for a single, airfoil, which can be 
solved by known methods. 

The analysis involves the assumption that the gap between a.ir~ 
foils is large compared to the chord length. In particular expres- 
sions, relating the characteristics of a cascade airfoil to those of 
a free airfoil are expanded in a power series in c/g, where c is 
the chord and g the perpendicular distance between airfoil chord 
lines in the cascade, and powers of c/g higher than the second are 
neglected. Definitions of the symbols used are found in Appendix A, 


THEOBY 

Consider an infinite unstaggered cascade of identical two- 
dimensional airfoils, as represented in figure 1. The cor.figura.tion 
is specified by the chord/gap ratio c/g, the camber— line shape and 
thickness distribution of the individual airfoils. For an indivi- 
dual airfoil, the incident velocity is V with a corresponding 
density of p and angle of attack a' . 

The pressure distribution on a typical cascade airfoil will 
be compared with that which would be obtained over a single airfoil 
of the same shape in a stream of velocity V and density p. The 
angle of attack a of the single airfoil will not, in general, be 
the same as the angle of attack a’ of the airfoil in cascade, and 
an expression relating the two angles will be given. In the analysis 
to follow, primes will be used to designate properties of the cascade 
airfoil. 


Aerodynamic Characteristics of a Cascade 
Airfoil of Given Profile Shape 

Effect of camber. - The chordwise lift distribution for the free 
airfoil is given by the Kutta-J oukowski relation 
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where L is lift per unit span, r is the vorticity per unit 
span, and x is the chordwise distance from the leading edge to 
the point in question. Ihe lift-distribution coefficient P is 
then defined as 


„ 1 dL 2 dr 

P - q 3x s 7 lx 


( 2 ) 


In Appendix B it is shown that the stream velocity at the 
cascade airfoil is uniform so that the lift distribution for the 
cascade airfoil may be expressed as 


dL* _ dr' 

sr = pY sr 


(3) 


and the lift-distribution coefficient of the cascade airfoil referred 
to the dynamic pressure q is 


1 dL' _ 2 dr 1 
q dx ¥ dx 


w 


From equations (2) and (4), 


p* _ p - 2 

V 


\ dx dx / 


(5) 


It is now desirable to express equation (5) in a more convenient 
form. It is assumed that the vorticity distributions of the free 
and cascade airfoils, respectively, may be represented by the 
following series (references 1 and 2): 


££ = 2Y 
dx 



cot ie + 


dr 1 

dx 


2Y 



cot ^0 


A n sin nB 


n~i 


+ V A n f V-9 ) 

L, ‘ J 

n- 1 ' 
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The new parameter 0 is related to x by the equation 


x = & (1 - cos 0) 


(7) 


Equation (5) then becomes 


CO 

P* - P » k [ (A 0 ‘ - Aq) cot g0 + y (A n r - A n ) sin n0 ! 


( 8 ) 


n=r 


The coefficients in equation (8) will now be evaluated by 
considering the conditions of flow at the airfoil boundaries. Let 
v and v* be the vertical components of velocity induced, respec- 
tively, by the free airfoil and a. particular cascade airfoil, and 
let Av be the vertical component of velocity induced by the other 
airfoils of the cascade. For small angles of attack, in order that 
the flow be tangent at the surfaces, the slope dy; /dx of the air- 
foil camber line (which is the same for both free and cascade air- 
foils) must be given by the following relations: 


dye 

dx 


V 

a + T 


dye 

dx 


= a* 




(9) 


Appendix B gives the vertical components of velocity induced 
by the vorticity and source— sink distributions which represent the 
airfoils of a cascade. For simplicity, the following symbols are 
used : 


a 


_ Tt 2 C 2 


(Bll) 


X » Jl-U 3 


(B4) 


where M is the Mach number of the undisturbed stream. 



NACA EM No . A7E29 


5 


The velocity components calculated, in Appendix B are 

2 = \ A 0 + ^ A n cos n 0 

n=i 
00 


y — 

”• — \ I — Ao * + Ap * c os n0 

v V u 


(B5) 

(B 6) 


n-i 


and neglecting terms involving a to the second and higher powers 
f 1 = ~2 SL [(Ao* + ^Ag*) - (2 Aq* + Ax’) cos 0 ] (B25) 


Substituting these relations in equations (9), there is 
obtained 


&yc _ 

dx “ 


d 7c 

dx 


V 


-Ao + 

^A n cos 

n0^1 

J 




n=i 




(" “Ac 1 

00 

+ Va • 

+ ) K n 

cos n© J 1 

i 

(10) 

\ 

L — j 
n= i 

/ ! 




-2 J [ (A 0 * + $Ag‘) - (2A 0 * + A x «) cos 0 ] 


J 


Since equations (10) are equal trigonometric series, the 
coefficients of cos n0 can *be equated: 

a — AAo = cc T — AAo* — 2 ^ (Ao* + "gAg * ) * 

a ' i 


\Ai m XAi« + 2 a (2A 0 r + Ax') 

A 

Ag ss A^ 1 


V (u) 

f 


An = An’ (n £ l) 
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From the first of equations (ll), it is seen that the quantity 
(Aq 1 - A 0 ) depends on the angle of attack a of the free airfoil. 
This angle may he arbitrarily defined so that 


A 0 * - A 0 = C 


( 12 ) 


Combining equations (11) and (12), the following relations are 
obtained: 


A 0 « — Ao — 0 

Ax' - Ax = -2 (SA 0 * + Ax') 

A 2 * - A s - 0 (13) 

An* — An ** 0 (n £ l) 


and 


a* — a = 2 ^ (Ao* + \As ' ) (1*0 

Making use of equation (13), equation (1*0 may also be written as 
follows: 

a’ — a, = 2 j- (Aq + ^As) (15) 

From equations (8) and (13), there is obtained 

p* _ p - _8 JL (2A 0 ' + Ax') sin 9 (l6) 


Changing to the unprimed coefficients with the use of equations (13), 
and neglecting terms in a s } the expression becomes 


p* — P *= -8 At (2Ao + Ax) sin 9 


(IT) 
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It should not be Interpreted from equations (l6) and (17) 
that Ai* and Ai are equal. Equation (17) is the result of 
dropping terms in 0 s and not of equating Ax’ and Ax. The 
difference between Ax' and Ax to the order a is still given 
by equation ( 13 ). 

The section lift coefficient for the free airfoil is 

Cl= i' lpd (i) (18) 


Substituting from equations (2), ( 6 ), and (7) and performing the 
integration, one obtains 


C 7 ** « (2A 0 + Ax) 

Similarly, the quarter-chord moment coefficient is 

Cm c/4 = P (k ~ f ) d (f) = “ | - A s) 

From equations (19) and (20) it follows that 

% 

2Ao + Ai = — “ 

X 3t 

Ao + sA 2 = (ci + 4 C m c / A ) I 


(19) 


( 20 ) 


( 21 ) 


With use of equations (21), equations (15) and (17) then become 


\ 

a' - cl = (ci + b c ffic/4 ) 


( 22 ) 
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and 



Effect of thickness .— From Appendix B, the horizontal velocity 
at an airfoil in cascade is greats?.’ than that at a. single airfoil 
"by the amount Au, which is given by the .relation 


An Ac 
V X 3 


(B24) 


The quantity A is a function of airfoil thickness given by 
equation (B20) and tabulated for various airfoils in table I, which 
has been taken from reference 3 . 

The increment of velocity given by equation (B24) can be 
added to the velocity of the undisturbed stream Yx to give the 
true incident velocity V: 


T = Yx + Au 

From equations (B24) and (25), 

Jx - i _ Ao 

Y X 3 


(25) 

( 26 ) 


The density pi far ahead of the cascade , is given by the 
following equation, derived for isentropic flow: 



where 7 is the ratio of specific heat at constant pressure to that 
at constant volume. Using the binomial expansion and neglecting 
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terms containing powers of [(Vi/V) 2 — 1] higher than the first, 
there is obtained to the order of a 


fii. - 1 + MB 
P 


( 2 ?) 


The ratio of the dynamic pressure q at the center of 
pressure to the dynamic pressure q^ far ahead of the airfoil 
is then found to be 


SL * 1 + 2 4 & 

qi X 


( 28 ) 


where 


i - l-£ 


(29) 


The lift-distribution coefficient P' referred to the 
dynamic pressure qi is given by 


- SL_ p* 

r qi ‘ 


(30) 


The cascade section lift coefficient is 


c x' = f P* d(x/c) = §- / P“ d(x/c) 
Jo yi J° 


(3D 


Substituting equations (24) and (28) in equation ( 31) and 
performing the integration, there is obtained, neglecting terms in a 2 , 

o\ x — ci = —2 ^ ^1 - (32) 

Equation (32) gives the relation between the lift coefficient of a 
cascade airfoil at angle of attack <%' and that of a single airfoil 
at angle of attack a. The relation between a* and a is given 
by equation (22) . 
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The Mach number M’ must now he related to the Mach number 
M for the corresponding single airfoil. The velocity of sound ai 
is related to a by 


ai _ l a 

a l *** f rn a 
1 


(33) 


For isentropic flow 


Ti = T -T 1 — M a 




which becomes., using equation (26) 


Ti = T 


i + (r - i) m 2 M 

A 


(34) 


Hence equation (33) becomes, to the order of a 


ai 


(35) 


and so 


M< -II. 
ai 



a 


a. x 


Using equations (26) and (3i?) then 


M* 


Ac 

"F 




(36) 


where 


h = 1 


+ ( 2 i i ) M2 


(37) 
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Determination of pressure coefficient.— The following coeffi- 
cients are useful in expressing the pressure on an airfoil surface: 


Si 


H— Pi . 


pi 


pi-p 

q 


(38) 


where is the local static pressure on the surface of the airfoil 

and H, p, and q are, respectively, the total head, static pressu 
and dynamic pressure of the undisturbed stream. The variation of 
with Mach number, assuming that the ratio of specific heats is 1.4, 
is given (reference 4) by the equations 


H ** p + q (1 + tj) 

jyr 4 

to 


i + n =i + ^ + t^ + 


_M 6 
l600 "* 


(39) 


From these equations. 


Si = 1 + rj - Pi (40) 

A graph of rj as a function of Mach number is given in figure 2 
The subscript l in equations (38) and (40) will be replace by U, 

L, or f to denote, respectively, the upper and lower surfaces of 
an airfoil and the surface of its symmetrical base profile. 

According to reference 5, the velocity Vf along the ha.se 
profile of a single airfoil is given, for incompressible flow by 
the expression 




where the pressure coefficients aro referred to 


( 41 ) 


Pu - 


2?; Pt- 


PL-P . 



q. That is. 


ixW 
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and p is tho static pressure of the stream corresponding to q. 

The velocity Vf ' along the "base profile of the airfoil in a 
cascade in terms of Y is simply 


Yf 1 _ Vf 
Y " Y 


so that 


J l~Pf * « V 1 "^" 


m 


since cascade base— prof ile pressure coefficient Pf* is referred 
to the dynamic pressure q and to the static pressure p rather 
than to the corresponding quantities in the undisturbed stream. 
Analysis (reference 3) she vs that equation (42) holds also for 
compressible flow, to an approximation of the same order of magnitude 
as others already made . 


The upper- -and lower— surface pressure coefficients referred to 
q a.re from reference 5* 



) ( *a) 


Although equations (^3) were derived for incompressible flow, it 
can be shown that they are applicable in the case of flow in a . 
compressible stream. 

The coefficients Eg* and Pj* can be converted into the 
coefficients Pu* and Pg* ref erred, to qi (the dynamic, pressure 
in the undisturbed stream) and the corresponding static pressure by 
means of the following equations, which are based on equations (28) 
and (40) : 
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su* = l + ti - Hr* 

s L * - 1 + tj - Pl* 

su» - £ su* a su* (i 

sl' » sl* - s L * (i 

Pu’ = 1 + ri* — Su f 

pl* = i + V - sl' 


when tj * is the value corresponding 

Determination of Pressure Distribution for 
a Given Profile Shape in Cascade 

A method can now be outlined for finding the pressure distri- 
bution over an airfoil of given profile in an unstaggered cascade. 

It will be assumed at the outset "that the Mach number M* of the 
approaching air is known and that the desired lift coefficient c^ 1 , 
of each of the cascade airfoils is given. 

1 . The Mach number of the corresponding free airfoil, M, is 

determined from equation (36) to a first order in a as 



J 

M* given by equation (36). 


M = 


0 



( 45 ) 


wherein 



X = ft 1 ) 2 
p = 1 + i (M r ) 2 


(Bll) 

( 46 ) 


( 4 ?) 
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and A may "be obtained from table I or from equation (B20) . 

2. The lift coefficient of the corresponding free airfoil, 
ci, is found from equation (32) by neglecting second order terms 
in a as 



(48) 


wherein | is given by equation (29) . The angle— of-attack relation, 
in degrees, is 


a* = a + 57.3 (ci + 4 cm Cy 4) (22) 


3 . P = Pj, — Pg for compressible flow at the Mach number M 
is found for the single airfoil at the lift coefficient' c i. 
Preferably an experimental pressure distribution at the appropriate 
Mach number and at approximately the same Reynolds number should 
be used but a. theoretical distribution (e.g. , by references 2, 6, 

7, or 8) modified by the Glauert— Prandtl on Karmdn—Tsien rule is 
satisfactory (references 9 and 10) . 

4. P* is found from equation (24). 


5. a/1— Pf* is calculated from the free airfoil pressure 

distribution from equations (4l) and (42) 

J 1-Pf * =V /1-Pf = t (49) 


6. (1— Pf*) and P* are combined to give ( 1— Py* ) and 
(l- Pi*) from equations (43) 
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1? 


1 _ p..* „ ♦ fr) 

i _ P * = (i-Pf» -£p*) ; 

i-Pf* 


2 


m 


7. The cascade pressure coefficients Pg* and Pj/ are then 
found by means of equations (44): 

Sg* 

S L * 

Sg» 

Sl' 

Pg’ = 1 — Sg* + q 1 

. pl* * i - s L ' + n* > 

where rj a.nd q* corresponding to M and respectively are 

found from figure 2. 


1 - Pg* + q 

1 - Pl* + T] 

A erg 'N 
X 3 7 


= 1 + 2 


Su* 


(i + 2^| ) s L * 


) 


(44) 


Determination of Profile Shape From a. Given 
Pressure Distribution In Cascade 

The procedure for obtaining the profile shape for a. given 
pressure distribution in cascade involves a. process which is essen- 
tially the reverse of that Just outlined. First the required pres- 
sure distribution is drawn. The initial choice must be made skill- 
fully with reliance on experimental pressure distributions so that 
the distribution chosen shall correspond or nearly correspond to 
that obtainable with a real airfoil. From the chosen distribution 
the coefficients cj‘ and cmcy^’ are determined by graphical or 
numerical integration and the pressure coefficients Pg* and Pl* 
afo road off at selected points. The method then consists of the 
following steps: 
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1. Pu* and Pl* are found by applying equations (44) in • 
reverse order, neglecting terms involving a to the second and 
higher powers. 

Su* * 1 + V - Pu* -s 

Sl* » 1 + V - P L » 

su* =f 1 - 2^1^) SuK 

. ; " ■ > (50) 

SL* = (l - 2^ j SL* 

?U* = 1 + rj — Sy* 

pl* = i + t) — sl* > 

2. (l-Pf) is found by means of the following equations, which 
are readily obtained from equations (49) 


Vi-Pf = V "i-Pf* 


*/ l—Pij* + ^ 


(49) 


The base profile should now be checked to see if it satisfies the 
closing condition given in reference 2. If the assumed pressure 
distribution does not correspond to a closed shape, it must be 
slightly modified until it does . 

3. The single-airfoil lift-distribution coefficient P is found 
from the relations 



4. The final upper— and lower— pressure coefficients are calcu- 
lated by means of the following equations from equation (43): 
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IT 




1 - 


(l-Pf +?P) S 


PL = 1 - 


l~Pf 


(1_ P -ip) 2 


1-P* 


( 53 ) 


5. Tlie angle of attack of the single airfoil is then given in 
degrees by 


a = a' ~ 57-3 ^ (c* 1 + 40^^) 


( 54 ) 


which is a modification of equation (22). 

The problem of finding a profile shape which will have a given 
pressure distribution in cascade has now been induced to the analogous 
problem for a single airfoil. It is now possible to correct the 
pressure distribution given by equation (53) by the method of 

reference 10, or simply by multiplying Pu and Pl by */ 1-M 2 . 

An airfoil shape can be designed to give this pressure distribution 
by the method of reference 2. In the application of this method, 
it is generally necessary to make some changes in the pressure 
distribution, as previously noted, in order to make it correspond to 
a possible distribution for an actual airfoil. If the initial distri- 
bution has been well chosen, however, the changes will be minor. 


The Choking of Cascades 

In the compressible adiabatic flow of a fluid in an elementary 
stream tube of varying area ki, the mass flow must be constant so 
that the logarithmic derivative 

d-P; + 

PI Vi A* 

must vanish. Now the density pi and the velocity Tj are related 
to the pressure p^ by Bernoulli’s equation 
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&PI _ <3-P l x dp i 
Pi d P l Pi 




The quantity dp /dp ^ is, of course, the square of the velocity of 
sound so that Bernoulli's equation my he written 


dp 2 Vim 

Pi a.2 2 


= -Ml 


2 

V l 


where is the local Mach number. Combining this expression 

with that of the logarithmic derivative then 


( 1 -^ 2 S ) 


m dfi.2 

Vl Al 


From this relation it is seen that for a subsonic flow the 
area must decrease for a velocity increase, while at supersonic 
speeds the area must increase for a velocity increase. When the 
Mach number is unity, then dA = 0, so that sound speed is only 
attained where the area is a minimum. 

Considering the flow through the ca,scade of figure 1 as 
essentially unidimensional, then it is apparent that if the flow 
past the plane AA attains sonic speed, the mass flow through the 
cascade cannot be further increased and the cascade flow may be 
said to be "choked." Of course, the flow through a cascade is not 
unidimensional, but experience with the similar phenomena of choking 
in wind tunnels has indicated that the assumption of unidimensionality 
of flow yields calculated choking Mach numbers in good agreement 
with experiment. 

For a unidimens ional flow it is shown in reference 3 that the 
ratio of the free area of the undisturbed stream A to the minimum 
flow area A m is related to the choking Mach number of the free 
stream M c h by 
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For the approaching stream which will pass between any two airfoils 
of tho cascades the free cross-sectional area is 


A = g cos a' 


while the minimum area between the two airfoils is 


Am 



and hence 


JL 

Am 


- cos a* 
T t 'o 

-A — 77 


(56) 


where t is the maximum thickness of an airfoil. 

Properly the thickness of the boundary layers on each surface 
of the airfoil should be added to the geometric thickness to obtain 
an effective thickness t e . Using this value with equations (55) 
and ( 56 ) and a value of 7 (for air) of 1.4 then 



J'4 c h C 03 a* 

V' 5 ) 


(57) 


A practical estimate of the choking Mach number may be found by 
assuming the angle of attack to b© so small that cos ex’ = 1 . and 
that the effective thickness Is the geometric thickness . With these 
assumptions , values of the thickness gap ratio as a function of the 
choking Mach number are given in figure 3 . 

It is shown in reference 3 that there exists another possibility 
that choking may occur in the wake of the airfoils for very thin 
profiles as the result of the action of viscosity. In most cases 
this type of choking will not be of practical importance . 


DISCUSSION 

In order to check the accuracy of the equations which have been 
developed, pressure distributions for incompressible flow have been 
calculated at three lift coefficients for an unstaggered cascade of 
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NA.CA bkl2 airfoils. These cases were chosen to permit comparison 
with those which have been determined in reference 11 by the method 
of conformal transf ormation. The chord/gap ratio in all cases is 
1.03 and the lift coefficients considered are 0, 0.5* and 1.0. This 
comparison subjects the approximate theory of this report to a 
rather severe test since the analysis has been developed on the basis 
of a chord/gap ratio which is less than unity. The comparison with 
reference 11 indicates the agreement is good as may be seen in 
figures 4, 5* and 6. The single— airfoil pressure distributions, as 
obtained by conformal transformation in reference 11, are shown on 
the figures for the same lift coefficients as for the cascade air- 
foils. (In Appendix C the calculations to obtain the pressure 
distribution of figure 6 are given in detail. Table II gives all 
the necessary computations and serves to demonstrate the simplicity 
of the method.) It is evident from the figures that one effect of 
cascading airfoils is to impose a "negative camber influence" on the 
pressure distributions as the lift coefficient is increased. Accord- 
ingly, for airfoils for use in cascade the camber must be exaggerated 
if a certain desired camber effect on the pressure distribution is 
to be obtained. It is of interest to note that the ca.lcula.ted angles 
of attach are also in reasonably good agreement as may be seen in the 
f oil owing tabulati on : 


Airfoil in Cascade 


Airfoil Alone 


£il 

a* (from 
reference 11 ) 

a' (method 
of report) 

a 

0 

_<=; o° 

—6*0° 

X 3 1 

• 5 

+1.8° 

+2.3° 

- 0.1 

1.0 

+9.7° 

+10.4° 

•4*4 , 0 


The corresponding ratio of the mean lift- curve slope of the 
airfoil in cascade to that for the airfoil alone. 


(&Ci/ /da* ) ca , 3 Ca ,de 
(&CL/da)i one airfoil 


is 0,53 by the method of reference 11 and 0.6l by the method of this 
report. 
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Unfortunately no pressure distributions over unstaggered 
cascades at high airspeeds could be found so that the validity of 
the compressibility corrections developed in this report could not 
be determined. They have been developed, however, in much the 
same way as the compressibility corrections of reference 3 * which 
havo since been found to be in good agreement with experiment. It 
is clear, however, that a3 the Mach number is increased the accuracy 
of the calculations will diminish unless the chord/gap ratio is 
simultaneously decreased. This will bo particularly true as the 
choking Mach number is approached. 

One matter of interest concerns the effect of compressibility 
on the lift— curve slope of a, cascade of NACA 4412 airfoils having 
cbord/gap ratios of 0, (i e., lone airfoil) 0,5, and 1.0 which is 
shown in figure 7. It is soon that the importance of the interfer- 
ence effects of the airfoils in cascade increase so rapidly with 
Mach number that, contrary to the usual expectation, a decrease in 
lift— curve slope with Mach number is indicated for high solidities. 

An examination of the equations for the choking of cascades 
indicates that care must be exercised with an unstaggered cascade 
to keep the airfoil thickness small for high solidities if relatively 
high free— stream Mach numbers are employed. Even for the cascade 
considered in figures 4, 5, and 6, for example, the choking Mach 
number as obtained from figure 3 is only about 0.64 based on the 
geometric thickness of the airfoils. 

The existence of a. boundary layer on the airfoil surfaces 
would, of course, increase the effective thickness which would 
reduce the choking Mach number. When the pressure gradients are 
strongly adverse, as in the case of the airfoil of figure 6, the 
boundary— layer growth will be increased and its effect on the 
choking Mach number will be more pronounced. The proper choice of 
camber consistent with the design lift will serve to reduce the 
sharp pressure peak, thereby improving the critical as well as the 
choking Mach number. 


Amos Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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APPENDIX A 
'List of Symbols 

The following symbols are used through this report: 
c airfoil chord 

g distance perpendicular to chord between airfoils 

in cascade 

A factor depending upon shape of base profile (See 

equation (B20) and table I.) 

o factor depending upon solidity of cascade jjjg- ~ 

a angle of attack, degrees 

c i section lift coefficient 

Cja c / 4 section quarter— chord-moment coefficient 

V velocity 

M Mach number 


choking Mach number 

compressibility correction factor^*/ 1-M 2 j 
compressibility correction factor ^1 — 

compressibility correction factor ^ 1 + M 2 ^ 

ratio of specific heat at constant pressure to 
specific heat at constant volume (c- 0 /cv) 

total head 


p static pressure 

q dynamic pressure 


P 


mass density 
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n 

p 

?Z 

Sz 

x 


r 

<f> 


T 


yt 


t 

dy c /dx 

0 , 

ar/dx 

u 

V 

Av 


B n 

r 

L 


compressibility factor (See equation (39) and fig. 2.) 

local lift at any chord station in coefficient form 

local pressure coefficient (See equation 3&) 

local pressure coefficient (See equation 38 .) 

coordinate of points on chord line as measured from 
leading edge 

angular coordinate of points on chord line (See 
equation ?.) 

radial distance in polar coordinates 

polar angle in polar coordinates (positive counter- 
clockwise) 

inclination of the radius r to a line normal to 
the airfoil chord (See fig. 1.) 

ordinate of base profile 

maximum airfoil thickness 

slope of mean— camber line 

source strength 

vorticity per unit length of chord line 

horizontal component of velocity 

vertical component of velocity 

total induced vertical velocity at the airfoil 
under consideration due to the other airfoils 

Fourier coefficients (See equation 7.) 

Fourier coefficients (See equation (331*0.) 

circulation per unit span 

lift per unit span 
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T absolute temperature 

Superscripts 

1 used, to distinguish properties of cascade airfoil 

* donotes cascade airfoil characteristics as coefficients 

referred to dynamic pressure q of incident stream 
at center of pressure 


Subscripts 

1 denotes values in stream far ahead of cascade 

l denotes local conditions at point in fluid 

L denotes values on lower surface of airfoil 

U denotes values on upper surface of airfoil 

refers to "base profile 


f 
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APPENDIX B 

VELOCITY COMPONENTS INDUCED BY AIEFOILS IN CASCADE 


For inc empress ib le flow over a single airfoil, the vertical 
velocity (perpendicular to the chord) induced at the point on the 
chord x 0 hy its own vorticity distribution is 


v 


j r c M 

2nJ 0 x-x Q 


(Bl) 


To obtain the velocity for compressible ’ flow, the factor 


a/ 1-M g 
1-M 2 sin 2 $ 


(B2) 


must be applied to the right-hand side of equation (Bl) as is shown 
in reference 9* Here M is the Mach number of the undisturbed 
stream and § is the angle between the stream direction and a. 
line drawn from the vortex to the point in question. For a single 
airfoil $ is always close to 0° or l80° and the induced velocity 
in a compressible stream is approximately 

/r~~2 P C SjfE dx P° &£ dx 

v s Vl-Mf/ dx = 2u dx ( B3 ) 

2tc Jo x— x 0 2nJo x— x D 


where 


X 


m 


Equation (B3) can be integrated by substituting for df/dx and 
for x the following: 


4E 

dx 


= 2V 


CO 

^A 0 cot ■§£ + ^ A n sin ne^ 


n=i 


( 6 ) 


X S | (l - COS 0) 


( 7 ) 
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The details of the integration are given in reference 2, page 4. 
The resulting expression is 


v, 

V 


x(-A 0 





cos n0 


n= l 


(B5) 


The vertical velocity at a point on the chord of an airfoil 
in cascade is made up of three parts. These are (l) the velocity 
v* induced by its own vortices/ (2) the velocity Av v induced by 
the vortices of the remaining airfoils of the cascade, and (3) the 
velocity Av a induced by the sources and sinks of the remaining 
airfoils . 


By analogy with equation (35) : 




CO 





cos 


U6y' 


(B6) 


The vertical velocity induced by the vorticity distribution of the 
other airfoils is given by the expression 


AVy 


iV 

m-i 


AC 


aT* 

dz 


sin t 


m 


Jr 




( 


JTW 

1-M 2 COSS V 


m 


) 


4x 


(B7) 


The s umma tion is over the remaining airfoils of the cascade. The 
angle t is measured clockwise between the vertical and a. line from 
the vortex to the point in. question, as indicated in figure 1. The 
factor in parentheses is that necessary to give the correct result 
for compressible flow. It is derived from the factor (B2) with 
use of the relation sin ® = —cos t. Since rm = (mg/cos T m ) and 
tan v m =(x-xo)/mg so that 


and 


sin 


r m ~ 


X-Xp 

mg 




m 




COB T 
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the preceding equation may "be written 



Expanding the term in "braces in a power series in (x— Xo)/mg„ 
neglecting terms involving the cube and higher powers, and noting 
that 


00 



there is obtained 


Avv - 5 tv (x - xo) 4x (B9> 

If equations (6) and (?) are substituted into equation (B9) 
and the indicated integration is performed, there is finally- 
obtained 

~ 2 - [(A 0 * + |A 2 ') - (2A 0 ' + Ai«) cos 6 ] (BIO) 


where 


d-it§ < bii > 

It is Interesting to note that equation (BIO) can be more 
directly obtained by replacing each of the airfoils of the cascade 
by a single vortex. The important point in such a derivation is 
the proper chordwise location of this vortex at the airfoil center 
of pressure. 

In a. manner similar to that in which equation (BIO) was 
obtained it is possible to show that neglecting terms involving 
fc— Xc)/mg to the third and higher powers. 


= 0 

V 


(B12) 
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In appendix B of reference 3 the velocity components induced 
"by a source at the origin are given as follows: Hie velocity 

component perpendicular to the stream Av s is 


Q a/i-M 2 (sin cj>) 

2nr L 1-M 2 sin 2 5> 

and the induced velocity component parallel to the stream Au s is 

cos $ 

\/l~M 2 (l~M 2 sin 2 <S>) , 


J9L| 

2«r 


where Q is the mass flow divided "by the free— stream density. 

The vertical velocity induced "by the source— sink distributions 
of the remaining airfoils of the cascade can then he expressed as 


HI— 4t» 



m~x 


and the corresponding horizontal induced velocity is* after noting 
that cos <p = — Bin 


Au.3 = 


m=+co 

iy 

”„4 


r° 

dO.' /dx 

sin T m 

Jo 

I'm 

(l-M 2 cos 2 tjb) _ 


dx 


whene —ll dx is the source strength over an element of chord dx 
dx 

of the cascade airfoil. Substituting for rjn and cos t m and 
sin Tm as before, expanding, and neglecting the third and higher 
powers of (x-ab)/:mg, there is obtained 


Ava = 0 


re 


(B13) 


* - Sfei af (x - Xo) “ 
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The source-sink distribution is now expressed by the following 

series: 

= 2V (Bo* cot sin n0) (Bl^) 

n-i 

On substitution of equation (Blk) the integral in equation (B13) 
ca.n be evaluated,, obtaining., analogously with equations (BIO) 
and (B12) : ' 

—a - 0 (B15) 

and 

”3/^ = [ (Bo* "^irBs 1 ) — (2Bo* + Bx*) cos 0] (Bl6) 

Since the airfoil must be a. closed figure/ it follows that 

n c 

/ dx = 0 (B17) 

Jo dx 

By substitution of equation (Bl4) in equation (B17), there is 
obtained 

2B 0 r + Bx 1 =0 (Bl8) 

so that 

« ~g (2B 0 * + B 2 ') 

From thin airfoil theory (reference 2) it can be shown that 

2B 0 « + B 2 « = & [ & d (j) (B19) 

where y-j. is the ordinate of the base profile as measured from the 
axis of symmetry in terms of the airfoil chord at the corresponding 
x station. 

In the problem of determining the wall interference in a two- 
dimensional wind tunnel, an evaluation of the influence of the 
image base profiles yields an identical result if the limitations 
of thin airfoil theory a.re presumed. Lock (discussed in reference 3) 
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has evaluated the image base- prof ile effects as regards wall inter- 
ference in the incompressible case for profiles not necessarily 
thin and has found that (in the notation of this report) 


A = 



( 1 -Pfl) | 1 


~Tdv 


£t f 1 a/ s) 
cGc/ J W 


(B20) 


is a more precise value to replace 2Bo' + B 2 S where Pfj i3 the 
base— profile pressure coefficient in incompressible flow. Values of 
A computed by equation (B20) for various airfoils are given in 
table I which was taken from reference 3- Accordingly., as in 
reference 3 ? this value is used in this report so that the velocity 
increment due to the effect of sources and sinks becomes 


Aug _ Ac 
V " 13 


(B21) 


It should be noted that properly Lock's result for incompressible 
flow is 


^ u s _ . 
Y- ~ Ac 

V X 


(B22) ’ 


However., the expression for V is 


so it follows that 


V = Vi + Au 


Vi 

V 



Thus fox* incompressible flow 


Aus 

V 



(B 23 ) 


which to the first order in c is still equation ,(B22) so that 
equation (B2l) for compressible flow follows. From equations (B12) 
and (B2l) it is clear that the influence a.t the position of one 
airfoil of the cascade due to all other airfoils of the cascade is 
simply that promoted by the source-sink system 
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£& - Ace /’•r oIl'i 

V X 3 {B } 


Moreover, because Av a (equation (B13)) is zero, it follows from 
equation (BIO) 


&C. . -2 3. [(A 0 ‘ + |-A 2 «) - (2A 0 ’ + A x «) cos S] (B 25 ) 

V A 
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APPENDIX C 

SAMPLE PFJEESUS^ CALCUIATION 

Table II gives all the necessary calculations to arrive at the 
cascade pressure distribution shorn in figure 6. It was desired to 
obtain the pressure distribution for a cascade of NACA 4412 airfoils 
given the values ci 1 « 1.0 and c/g =* 1.03 for comparison with the 
distribution calculated by the method of conformal transformation 
given In reference 11. 

Since the assumed flew is incompressible, then X, P* and | 
are unity and q and r] 1 are zero. From table I, A is 0.237. 

The parameter a is 0.220, The lift coefficient for the correspond- 
ing airfoil In free air Is obtained from equation (48) as 1-335. 

By potential theory for the NACA 4412 airfoil the free— air angle of 
attack for ci « 1.335 Is 6.8° and cm c / 4 is —0.11 so that for the 
airfoil in cascade, by equation (22), a* = 10.4", 

Values of the pressure coefficient for the isolated airfoil 
were calculated from the free— air values of reference 11 using the 
method of reference 5- From these pressures P was determined and 
the P* was obtained from equation (24), P* and 1— Pf* (equal 
to 1— Pf) were then combined by equations (43) to give 1— Pg* 
and 1— Pi*- The cascade pressure distribution was finally obtained 
using equations (44) . 
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TABLE I.- VALUES OF A FOR VARIOUS BASE PROFILES 
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TABLE II.- CALCULATION OF CASCADE PRESSURE DISTRIBUTION OF FIGURE 6 (SEE APPENDIX C) 

[o = 0.220; A= 0.237] 
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Figure J.~ Choking Mach number as a function of the 
ratio of airfoil thickness to cascade gap . 
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Figure 4 Pressure distribution for NACA 44i2 Qirfoi! 
a tone and in cascade for - L03 and c % ' ~0, 
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Figure 6. - Pressure distribution for NACA 4412 airfoi! 
alone and in cascade for -y- - /. 03 and c£ - i.O. 
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